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Q ■ Abstract 

J>-j' The preprocessing stage of Shor's algorithm generates a class of quantum states referred to as 

: 

periodic states, on which the quantum Fourier transform is applied. Such states also play an 
' important role in other quantum algorithms that rely on the quantum Fourier transform. Since 

(N 

entanglement is believed to be a necessary resource for quantum computational speedup, we analyze 
Q_i' the entanglement of periodic states and the way it is affected by the quantum Fourier transform. 

, To this end, we derive a formula that evaluates the Groverian entanglement measure for periodic 

a 

^ I states. Using this formula, we explain the surprising result that the Groverian entanglement of 

the periodic states built up during the preprocessing stage is only slightly affected by the quantum 
. Fourier transform. 
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I. INTRODUCTION 



Quantum algorithms offer a potential speedup over classical algorithms in solving a num- 
ber of problems. The origin of this speedup is not yet fully understood, but quantum 
entanglement is believed to play a crucial role l|-|3[. Therefore, it is of interest to analyze 
the entanglement of the quantum register during the operation of quantum algorithms such 
as Grover's search algorithm and Shor's factoring algorithm ^-[sj. Currently, all known 
quantum algorithms presumed to provide an exponential speedup over their classical coun- 
terparts rely on the quantum Fourier transform (QFT) The most notable among them is 



Shor's factoring algorithm jlOl. During the operation of these algorithms, the quantum 
states of the register are characterized by multipartite entanglement. Unlike the case of 



bipartite entanglement 
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13|, the multipartite entanglement in a register of g > 2 qubits. 



is not as well understood, partly because no analog of the Schmidt decomposition was found 
for multipartite systems. 

In order to evaluate the entanglement of the state of a quantum register, an entangle- 



ment measure is needed 



14j-|l7|. Axiomatic considerations have provided a set of properties 



that entanglement measures should satisfy 
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17| . These properties include the require- 



ment that any entanglement measure should vanish for product (or separable) states. It 
should be invariant under local unitary operations and should not increase as a result of any 
sequence of local operations complemented by only classical communication between the 
parties. Quantities that satisfy these properties are called entanglement monotones. These 
properties provide useful guidelines in the search for entanglement measures for multipartite 
quantum states. Entanglement measures based on metric properties of the Hilbert space 
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18| and on polynomial invariants 19|, l20j were proposed and shown to satisfy these 



requirements. Specific measures of multipartite entanglement include the average bipartite 
measure [21] the Groverian measure j22|, and the geometric measure I23I. A major difficulty 
in the evaluation of multipartite measures is that they involve a minimization of a compli- 
cated function in a high- dimensional space. As a result, there are no general closed-form 
expressions for these measures. 

n 

The Groverian entanglement generated in Shor's algorithm was analyzed in Ref. j7|. It 
was shown that the entanglement builds up during the preprocessing stage and that the QFT 
has little effect on the Groverian measure. This is somewhat surprising since, in general, the 
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OFT operator tends to generate highly entangled states when it is applied on product states 
7|. Furthermore, the superior efficiency of Shor's algorithm is attributed to the QFT, and 
since entanglement is considered a necessary resource for quantum computational speedup, 
one would expect that the QFT will induce it. It seems as though for the purpose of quantum 
speedup it suffices for the QFT to simply operate on a highly entangled register rather than 
generate entanglement by itself. 

The states generated by the preprocessing stage of Shor's algorithm are called periodic 
states. These states consist of an equal superposition of basis states whose indices take 
the form i = jr + I, where j = 0,1,2,..., r is the period and / is referred to as a shift. 
It was shown by numerical sirnulations that these states have the property of not being 
further entangled by the QFT [7]. In this article we explain this surprising property using 
an approximated formula for the Groverian entanglement measure of periodic states. 

The article is organized as follows. In Sec. II we present the Groverian measure. The 
periodic states generated by the preprocessing stage of Shor's algorithm are described in 
Sec. Ill and their entanglement is analyzed in Sec. IV. The effect of the QFT on their 
entanglement is considered in Sec. V. The results are discussed in Sec. VI and summarized 
in Sec. VII. 



II. THE GROVERIAN ENTANGLEMENT MEASURE 

The Groverian measure of a quantum state of q qubits is based on the maximal 
overlap that {tp) may have with any product state |(/?), with the same number of qubits. The 
smaller this overlap gets, the more entangled the quantum state becomes. We define the 
square of this overlap as 



PmaxW = , ^ , max , (1) 

where \ipm), m = 1, . . . ,q are single qubit states. This quantity cannot be decreased by local 
operations and classical communication between the parties holding the different qubits. 
Therefore, any nonincreasing function of Pmax that vanishes for product states (where Pmax = 
1) is a valid entanglement measure. Among all these possible measures, we have found it 



useful to use the logarithmic Groverian entanglement measure 



24| 
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G'(^) = -ln(P^axW), (2) 

to which we refer later in this article as the Groverian measure. This measure has three 
important advantages over other possible measures: (i) It is intrinsically a multipartite 
measure, rather than an average over bipartite measures for different partitions, (ii) It takes 
values in the range [0, oo), providing a better resolution than measures that are limited to 
the range [0, 1). This is particularly important in the case of highly entangled states with a 
large number of qubits. (iii) This measure is additive in the sense that if the subsystems A 
and B are not entangled with each other, then G{ipA ® ^b) = G{iPa) + G{iPb)- 

The problem with Groverian-type entanglement measures (and with many other proposed 
measures), is the difficulty involved in calculating them for general quantum states. The 
calculation of Pmax involves finding the product state for which Kyjl-?/;)!^ is maximal. The 
state \(p) is then referred to as the nearest product state. This is a maximization problem 
in a high-dimensional space. To evaluate the dimensionality of this space, we first consider 
the two-dimensional Hilbert space of a single qubit, which has four real parameters. 
The normalization and the insignificance of the global phase make it possible to express the 
quantum state of a single qubit in the form 



\iPm) = VT^^\0) + ^e''"^\l), (3) 

with only two real parameters. The first parameter, Xm, represents the balance between |0) 
and |1) in the corresponding qubit, and takes values in the range [0, 1]. When Xm = the 
qubit is in the |0) state and when Xm = 1 it is in the |1) state. The second parameter is the 
relative phase 6„i- Altogether, finding Pmax involves maximizing a suitable function in a 2q- 
dimensional space. This function typically exhibits a large number of local maxima. For a 
large number of qubits, this calculation requires significant computational resources, except 
for some special quantum states for which analytical formulas for Pmax niay be found. So far, 
a general formula is known only for two-qubit states (using the Schmidt decomposition) and 
br a very restricted set of highly symmetric states that generalize the GHZ and W states 
gI- Recently, some attempts were made to find formulas for the entanglement measure of 
various three-qubit states 25|, |26|. However, so far such formulas were found only for a 



restricted set of states. For arbitrary states of more than two qubits, the Groverian measure 
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can only be calculated numerically. 

A numerical scheme for calculating the Groverian entanglement is described in Ref. 
In each step of the scheme, one qubit 1 < mo < q is selected, and the parameter values for 
all the other qubits in the product state are fixed. The values of x^o and for which 
the overlap with lip) is maximized can then be found analytically. Repeating this step for 
every qubit several times, a maximum for the overlap over all the x^'s and 6'm's is found. 
Using such a formula to locate the maximum with respect to each qubit is much faster than 
successive evaluations of the overlap in a steepest descent method. This is due to the fact 
that the evaluation of the overlap requires resources that are exponential in q. Such a series 
of successive jumps in parameter space is also less likely to be misled by local maxima than 
the steepest descent method. A slight improvement to this scheme was used in Ref. j2^ . 
where at each step an analytical maximization was preformed over two qubits, using the 
Schmidt decomposition. This improved scheme is used in the present work as well. Still, 
the numerical calculation of the entanglement is time consuming. Furthermore, the lack of 
an analytical formula makes it difficult to achieve a better understanding of multipartite 
entanglement and its relation to quantum-computational speedup. Thus, it is worthwhile 
to search for analytical formulas for the Groverian entanglement of states that are relevant 
to quantum algorithms. Such an approximated formula is derived in Sec. IIVI 



III. PERIODIC STATES IN SHOR'S ALGORITHM 



Shor's algorithm aims to find a factor of a given nonprime integer A^. This is done by 



reducing the factorization problem to the order- finding problem [27|. In the order-finding 
problem one selects an integer y which is coprime to A^ and finds its order modulo A^, 
denoted r. By recalling that the order of y modulo A^ is the smallest integer such that 
y^ = l(modA^), one can see that when exponentiating y"'{modN) for a = 0,1,2,..., the 
resulting series will be periodic, with a period r. This can be done simultaneously for all 
values of a by constructing a superposition of the Q = 2^^ computational basis states in a 
quantum register with q qubits 
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The proper choice of q is described in Ref. 10|. The result of the modular exponentiation 
can be held in an auxiliary register: 



. Q-i 

Measuring the auxiliary register will randomly select one of its values, z = y\modN) for 
some < / < r, and will also filter out from the main register only those values of a for 
which = z. Since the series y'^{modN) is periodic in a with a period r, the values of a 
that remain will make out an arithmetic progression with a common difference r, and an 
initial term /. The main register will then be in a state we refer to as the periodic state of q 



qubits, with period r and shift I (following Ref. 28|): 

rn - n 

(6) 



j=o 

This ends the preprocessing stage of Shor's algorithm, and here the QFT is applied. In anal- 
ogy to the discrete Fourier transform (DFT), the QFT is used in order to reveal periodicities 
in its input [1]. In particular, the amplitudes of the state lip^i) ^^^ke out a periodic series, 
and when the DFT is applied to it, the resulting series can be approximated by a periodic 
series of the same sort, that is, one in which the indices of the nonzero terms make out an 
arithmetic progression. In the resulting series, though, the common difference is Q/r, the 
initial term is zero, and additional phases are added. This can be seen through the exact 
formula for the resulting series {yj)f^Q, given by 

^ 1 sin(7rjrA/g) .2.,[,+ i,(^„i)] 
VQZ sin(7rjr/g) ' 

Since applying the QFT to a quantum state is equivalent to applying the DFT to its ampli- 
tudes, the action of the QFT on periodic states can be approximately described as: 



^ I^^Mo)^ (8) 

where relative phases are ignored. Within this approximation, the QFT induces two changes 
in the periodic state, in analogy with the DFT: the period is changed from r to Q/r, and the 
shift is changed from / to (Fig. [T]). This removal of the shift is the crucial effect that makes 
it possible to extract the period in the next step, in which a measurement is performed. The 
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result of the measurement is not affected by relative phases; thus, they can be ignored. The 
measurement result is an integer close to jQ/r for some j, and dividing by Q we are left 
with a number close to j/r. A continued fraction expansion can then reveal j and r. If we 
had measured the state before applying the QFT, the result would have been an integer of 
the form l+jr, which does not allow finding the period without knowledge of the shift. It is 
clear that periodic states are central to the success of Shor's algorithm, and their importance 
is further demonstrated by their use in several other algorithms that make use of the QFT, 



all of which belong to a class of problems derived from the hidden subgroup problem 



28|. 



IV. THE ENTANGLEMENT OF PERIODIC STATES 



A. Equal superposition states 

The lack of a general analytical formula for the Groverian entanglement makes it hard to 
construct a model that explains the fact that the QFT does not seem to affect the entan- 
glement of periodic states. Nevertheless, an approximated formula may suffice, provided it 
remains close to the exact value when the number of qubits increases. Such an approximation 
may be attainable, since the set of periodic states is only a restricted set. 



1. Definition of the equal superposition states 

Let us first consider a somewhat less restricted set of states, namely, states that are 
superpositions of any number of computational basis states, with amplitudes that have 
equal magnitudes and zero phases. Given some nonempty subset of basis states S, we refer 
to their superposition as an equal superposition state of q qubits (ES state), 

l^^) = ^El^)' 0C^C{O,l,...,Q-l}, (9) 

V 1*^1 k&S 

where l^l is the size of the set 5*. Clearly all periodic states are ES states. Some ES states 
are nonentangled, like the computational basis states themselves. Another example is the 
complete ES state jr^), which is the superposition of all basis states and can be written as 
1+)®'?, where 

l+) = ^(|0) + |l)). (10) 



Other ES states are maximally entangled, such as the GHZ and W states j29|. 

Given any ES state, we would like to refer to the binary representation of each G 5* as 
k = ji . . .jg, where each is either or 1, and m = 1, . . . , g is the index of the qubit. This 
notation in hand, along with the notation of Eq. we can write down the overlap of the 
ES state with the product state \(f) as a function of the x^'s and Om's. It is this overlap 
function we need to maximize in order to find Pmax^ 



' ' ji...jqeS m=l 



where 



, \/T^^a^ 7 = , ^ 

C]" = <; (12) 

\/ j !• 

To be more precise, we need to maximize P = which is equivalent to maximizing the 

magnitude of the complex function /|, ignoring its phase. In fact, we can fix the relative 
phases 6^ to zero, since Eq. (ITT]) then becomes: 

V I I jl...jq&S m=l 

which is clearly not smaller in magnitude. This reduces the dimension of the search space 
from 2q to q (this is actually a special case of an observation already made in Ref. jol). 
Another simplification arises in an ES state for which there is one qubit, 1 < m < g, that 
all the basis states in the superposition "agree" on (that is, j.^ is the same for all k E S). 
This qubit is not entangled with the rest of the qubits and can be factored out. We shall 
call such an ES state reducible, since we can fix Xm = jm and reduce the dimension of the 
search space by one. 



2. Examples: Special high symmetry states 

We have so far reduced our problem to finding the values of Xi, . . . ,Xq for the nearest 
product state \ip), to a g-qubit nonreducible ES state {ipg), by maximizing the overlap 
function /| in Eq. (fT3|) . Let us consider two special cases of this problem, which were 
already analyzed in Ref. (6|. The first one is the g-qubit GHZ state: 
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For the GHZ state, the nearest product state is any one of the two basis states that comprise 
it: 



|(/p) = |0...0) or = (15) 

and for both states: Pjnax(GHZ) = 1/2. The second special case is the 2n-qubit balanced 
generalized W state, denoted by \(f){n,2n)). It consists of all the basis states of 2n qubits 
that have n zeros and n ones: 

_ 1 

\<f>{n,2n)) = (^^) ' Yl (16) 
In this case the nearest product state is the complete ES state 



\V) = \V), (17) 

and Pmax is given by 

Pma.[0(n, 2n)] = 2-2- ^ (18) 

3. The approximated formula 

Given a general ES state IV'I), our aim is to find the nearest product states jc/?). The 
results presented previously motivate us to examine two types of product states as candidate 
states: 

(i) \(p) = \k) for some k G S, which means that each Xm is either or 1. This yields 
P| = 1/1 5*1, where P| is the estimated value of P. 

(ii) \if) = \r)) (the complete ES state), which means that Xm — 1/2 for all m's. This 
yields PI = \S\/Q. 

Clearly, the first guess is better for a small 5", and the second guess is better for a large 
S. They become equally good for \S\ — y/Q. Thus, we can combine them into one improved 
guess: 



9 




\s\ 

M 
Q 



1 



s\<VQ 

S\>y/Q 



(19) 



We note that Pj is a lower bound on Pmaxi'^Ps)- We can also present a crude argument to 
support the claim that P| is a good approximation for Pmax('^s')- For a general product 
state \(p), consider expanding it to a superposition of basis states. In this expansion, we 
would like to maximize the number of basis states \k) that have corresponding basis states 
in lipg) (namely, for which k E S). The more such basis states, the larger the overlap will 
be (ignoring, for now, the amplitudes of the states). In the case of a small S, a single basis 
state is a good guess for \(f), since trying to vary any of the x^'s away from the edges of their 
range will add a lot of basis states to the expansion, and most of them will not be members 
of S, consequently decreasing P. In the case of a large S, a product state is desired with a 
lot of basis states in its expansion, since there are a lot of members in 5*, and in that case 
the complete ES state is hard to beat. 

4- A counterexample 

It turns out that the approximated formula presented above is not valid for all the ES 

j 

states. A state that violates this formula was analyzed in Ref. (6i]- This is the g-qubit W 
state, which consists of all the basis states that have q — 1 zeros and a single 1: 



Since for the simple W state 15*1 = q, and q < ^fQ for a large enough q (practically g > 4 is 
sufficient), our guess yields one of the basis states as the nearest product state and Pj = 
However, it turns out that the real nearest product state is 



which yields Pmax(W^) = [(? — 1)/Q']'^ • This is not only different from Pj; it is also asymp- 
totically different, as 




(20) 




(21) 
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(22) 



Pmax(W^) 




What property of the simple W state makes it violate the validity of the approximation? 
Note that for each qubit, the GHZ and balanced generalized W states have an equal number 
of zeros and ones across all the basis states. The simple W state obviously does not have this 
property. To understand why this is important, let us reconsider the maximization problem 
of finding the nearest product state \ip). 

We can explore the possible product states by taking a small variation around some basis 
state {ko). For two binary numbers ki and k2, let us denote the Hamming distance between 
them, which is the number of bits they differ on, as d{ki, ^2). We can then divide the set S 
to disjoint subsets according to the Hamming distance from ko'- 



Without loss of generality we can fix fco = (this can be arranged by applying local NOT 
gates which do not affect the entanglement nor the Hamming distances). The Hamming 
distance d{k, k^) is then equal to the number of ones in fc, and a small variation around k^ 
means that the x^'s are small. The terms in Eq. (fT3|) can then be grouped according to the 
subsets of S: 



The nth term in this expansion has n multiplicands of the form ^JXm, so it is dominant in 
respect to the {n + l)th term. Through this expansion we see that the nearest product state 
is in the close surrounding of |A;o) only if S contains a lot of terms within a small Hamming 
distance from k^. In the case of the simple W state, all the terms in 5* have a Hamming 
distance of 1 from k^ = 0, thus maximizing the term n = 1 in the expansion. This shows 
that a large value of the function fg can be obtained in the proximity of the state |0), and 
indeed this is the case. 

The case of the balanced generalized W state \(j){n, 2n)) is different. In this case, for each 
basis state, there are basis states at a Hamming distance of 2, and one basis state at 



S = S0iz}... U^g 

Sm = {k & S : d{k, ko) 



m} . 



(23) 




(24) 
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the maximal Hamming distance of 2n. In analogy, the GHZ state has a maximal Hamming 
distance between its two basis states. In both states there is no specific choice of ko for which 
all the basis states are within a small Hamming distance from it. This observation suggests 
that in general, for ES states that include basis states with large Hamming distances from 
each other, taking a small variation around some basis state does not aid in finding the 
nearest product state. 



B. Application to periodic states 

Returning to periodic states, the periodic state \ip'^i) is an ES state with: 

5 = {/,/ + r,...,/ + (A-l)r}; A = \S\ = ^— ^ . (25) 

r 

The overlap function in Eq. (fT3l) for this special case will be denoted by 

fU^i,...,x,) = ^ J2 11 CZ- (26) 

jl---jq&Sm=l 

The approximated formula for Pmax, which is the square of the maximum of the preceding 
function, becomes 

Approximating A ^ Q/r, we reach a simpler formula for G^, which is an approximation of 
the logarithmic Groverian entanglement of periodic states G{ip1.i): 

, - In i r < 
G'',= { ' ^ ^ (28) 

-ln§ r>v^ 

Looking at as a function of r (Fig. [2]) we see that it consists of two branches: 

(i) An ascending branch, which starts at when r = 1, and rises to a maximum at 
r = ^/Q. This branch covers the states with small periods, which have a large number of 
basis states, and for which the complete ES state is the presumed nearest product state. 

(ii) A descending branch, which starts at the maximum and descents back to as r reaches 
Q. This branch covers the states with large periods, which have a small number of basis 
states, and for which a basis state is the presumed nearest product state. 
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We shall now show evidence to support the validity of this approximation. We first note 
that in the case of an even period, the least significant bit has the same value for all the 
basis states in the superposition, which makes the state reducible. When factoring out the 
last qubit, the rest of the qubits constitute an arithmetic progression themselves, making 
them a periodic state with half the period. The value of the last bit depends on the parity 
of the shift: 

l^l.,.) = l<l//2j)®l^mod2). (29) 

Therefore, the problem of finding the entanglement of a periodic state with an even period 
can be reduced to the same problem with a state that has one less qubit and half the period. 
We shall concern ourselves from now on only with states that have odd periods. 

One important consequence of r being odd is that r is coprime to 2™ for every m. This 
means that the values of the mth significant bit in all the basis states go through a 2'"-long 
cycle with an equal number of zeros and ones. If the number of cycles is whole then the 
total numbers of zeros and ones in the mth bit are equal as well. However, in most cases 2™ 
does not divide A, the cycle is truncated and the equality is only approximate. Furthermore, 
the values of the n least significant bits make out an arithmetic progression with common 
difference r in respect to addition modulo 2", which means they also go through all their 
possible values in a cyclic manner. Therefore, the values of the different bits are in general 
uncorrelated. This means that for any basis state |A;o), the number of states in 5* that have 
n bits in common with it reduces approximately by half when n is increased by 1. This is 
because for each qubit m, about half of all the states have the same value for m as |A;o), and 
the same is true for any subset that is determined by the values of n qubits. We conclude 
that for a periodic state there is no basis state \ko) that more than any other basis state, 
has basis states in 5* within a small Hamming distance from it. This again suggests that the 
approximated formula for Pmax we have presented is vahd for periodic states, as opposed to 
the simple W state. 

An important observation concerning Eq. (12 7p is that for the descending branch, where 
the nearest product state \ip) is taken to be some basis state A; G S*, it is at least a local 
maximum of the overlap function f^.^. To see this, note first that there is no basis state 
k' & S that differs from k on exactly one bit. If there was such a state, the difference k — k' 
would be a power of 2, but this difference for a periodic state must be a multiplicand of r, 
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which is odd. For each 1 < m < q, the value of Xm is at the edge of its range (either or 
1), and trying to vary its value will turn \(p) to a superposition of \k) and \k') which differ 
only on the mth bit. Since k' ^ S the overlap function will necessarily decrease, and as this 
is so for all m orthogonal directions, \k) is a local maximum. The only question left then, 
for the descending branch, is whether \k) is a global maximum as well. 

Another justification to the approximated formula can be presented in the form of an 
induction, using the following recursive decomposition of periodic states. Considering the 
basis states that make up a periodic state, we can divide them into two subsets, according 
to the value of the most significant bit. Looking at the q — 1 remaining bits, we see that 
each subset makes up a component periodic state with the same period: 



Here, I' = —2''^^ (mod r) is the shift of the second component periodic state, Aq and Ai are 
the number of basis states in the component states I'lprl'^) I'^r'i'^) respectively. Clearly, 
the condition A = Aq + Ai is satisfied. When A is even, Ai = Ao = A/2. When A is odd 
Ai = Aq — 1, namely Aq = {A + l)/2 and Ai = {A — l)/2. Let us now assume that the 
approximated formula is correct for g — 1 qubits. If the states are in the ascending branch, 
then the complete ES state is the approximated nearest product state of the component 
states which make up the right-hand side of Eq. (130|) . Clearly this means that it is the 
nearest product state of the left-hand side as well, which is the periodic state of q qubits. 
On the other hand, if the states are in the descending branch, each of the component states 
has a (different) basis state as the approximated nearest product state. In this case, choosing 
one of these basis states gives the nearest product state to the state of q qubits. 

Finally, we present numerical evidence to support the approximated formula, described 
in Fig. [21 This figure shows the logarithmic Groverian entanglement of periodic states for 
some values of q, r, and /, computed numerically via the numerical scheme described in 
Sec. [TTl It also shows the entanglement according to the approximated formula: both the 
more accurate version given in Eq. (127]) and the less accurate version given in Eq. (1281) . As 
shown in the figure, the more accurate version of the formula agrees with the numerical 
results to a good precision on both branches. The less accurate version, however, does not 
follow the step-function-like behavior in the descending branch. This is clearly due to the 
formula for A, which includes a ceiling function that is smoothed out in the approximation 




(30) 
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\{Q-l)IA^Q/r. 



V. ENTANGLEMENT INDUCED BY THE QFT 

We have derived an approximated formula for the entanglement of periodic states with 
the hope of better understanding why the QFT does not increase their entanglement. First 
we note that this is indeed a special property of periodic states, since as illustrated in Fig. [3l 
the QFT operator in general changes the entanglement of quantum states. Looking again 
at Eq. ([H]), we can see that the QFT approximately takes a periodic state with period r and 
shift / to a periodic state with period Q/r and zero shift, up to relative phases. This was 
shown more rigorously in Eq. ([7]), where it was also shown that the relative phases previously 
ignored take a very special form. To explain this, we define a generalization of the regular 
ES state, by adding relative phases that depend on a parameter p: 

V 1*^1 kds 

In this phased ES state the relative phase of each basis state is proportional to the index 
of the state. An interesting property of this state is that it can be obtained from the 
corresponding regular ES state by local unitary operations, making it locally equivalent to 
the ES state: 

\^ e — j \Q e — j y e-Q^^' j 

Naturally, this means that both states have the same entanglement. As seen in Eq. (J?!), 
the relative phases added to the periodic state follow precisely the same pattern, with 
p = l+r{A—l)/2. Therefore, they can be ignored for all entanglement considerations. We are 
left with two changes the QFT induces in periodic states: removing the shift and changing 
the period. Interestingly, although the importance of the QFT in Shor's algorithm is in its 
canceling of the shift (which enables one to extract the period), this change is irrelevant to 
the state's entanglement, as the shift does not appear at all in the approximated formula. 
Therefore we are left only with the change of period. Looking again at the approximated 
formula for the entanglement of periodic states given in Eq. ( !28|) . we see that the values of 
the two branches, 1/r and r/Q, are swapped by the operation r — )■ Q/r. It is now clear why 
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= Gq^^ and therefore why G{ipri) ~ ^("^g/ro)- "^^^ QFT operator takes each periodic 
state in the ascending branch to a corresponding periodic state in the descending branch 
that has the same entanglement, and vice versa (recall that the QFT operator is its own 
inverse). 

Finally, we present numerical evidence to support the claim that the QFT does not change 
the entanglement of periodic states (for sufficiently large q). To this end we define the change 
of the Groverian entanglement of a state induced by the QFT 

AG{^) = G[QFT(^)] - G(^). (33) 

We examine this difference for periodic states as well as for random states. The random 
states are taken from a uniform distribution on the 2g-dimensional complex unit sphere. The 
Groverian measure of random states of q qubits exhibits a distribution that was calculated 
before for certain values of q Any unitary operator U (such as the QFT), when applied 
to a sample of these states, will produce states whose Groverian measure exhibits the same 
distribution. Looking at values of AG{iIj), for some states it is positive and for other states 
it is negative. The average of AGlip) over the random states is zero, but its distribution 
exhibits a certain width. In Fig. [3] we show the average of the absolute value of AG{iIj), 
which is an estimate of the width of the distribution of the values of AG{ip). For periodic 
states \AG{i/j)\ quickly decreases to zero with increasing q [Fig. |3]^b)], while for random 
states it changes only slightly [Fig. [3]^a)]. This demonstrates the special feature of periodic 
states, namely, that their Groverian measure is not affected by the QFT. 



VI. DISCUSSION 



The goal of the field of quantum algorithms is not only to find quantum algorithms 
that present a speedup over classical ones, but also to establish a deep understanding as 
to how this speedup is made possible. Currently, the best insight we can offer is that it 
is made possible by the combination of quantum superposition and quantum interference 
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3l| . Quantum superposition allows a sort of parallel computation, as all the states in 



the superposition go through the same unitary evolution in an independent manner. Since 
the results of each of the states cannot be accessed directly, we cannot fully exploit quantum 
parallelism. Nevertheless, through quantum interference the different parallel paths can 
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interact in a limited way, allowing us to access some global properties of the resulting states. 
Finding special situations where such a global property is the solution to some computational 
problem (like the period in Shor's algorithm is the solution to the factoring problem) is in 
fact the essence of quantum algorithm design. The role of entanglement in this model of 
quantum speedup is to allow quantum parallelism to reach its full extent, since product 
states cover a very small range compared with all possible superpositions 

To exemplify this model for specific algorithms, one must try to distinguish the role of 
quantum parallelism and quantum interference in each algorithm. In the case of Shor's 
algorithm, a superposition is built in the preprocessing stage, in a manner that does not 
make use of interference (modular exponentiation is performed on each computational basis 
state separately). This superposition in fact already contains the desired information (the 
period), and the QFT is merely needed to extract this information (by canceling the shift). 
Therefore, it can be argued that the preprocessing stage is where quantum parallelism is 
used, and the QFT introduces quantum interference. This claim is supported by the fact 
that the QFT does not increase the entanglement of the register. 

VII. SUMMARY 

We have shown that periodic states play an important role in Shor's factoring algorithm 
and pointed out that this is also true for other quantum algorithms that rely on the QFT. 
Focusing on entanglement as a necessary resource for quantum speedup, we set as our goal 
to explain the result presented in Ref. [7], according to which the QFT hardly affects 
the entanglement of periodic states. For this purpose we analyzed the entanglement of 
periodic states using the Groverian entanglement measure. We derived an approximated 
formula for the Groverian entanglement of periodic states and showed evidence to support 
it. Using this approximated formula, we presented a model that explains the aforementioned 
result. Finally, we argued that this result and the model that explains it strengthen our 
understanding as to the source of quantum computational speedup in Shor's algorithm and 
in general. 
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Figure 1: An example of the operation of the QFT on a periodic state. The periodic state of g = 8 
qubits with period r = 13 and shift / = 3 was transformed by the QFT. The amphtudes of the 
original state are given (a) , as well as the amplitudes of the resulting state (b) . The dotted vertical 
lines mark the multiplicands of the periods of each state (Q/r ^ 19.7). 
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Figure 2: (Color online) Entanglement of periodic states as a function of their period, r, for q = 8 
qubits and / = shift (a), q = 8 and / = 6 (b), q = 10 and / = (c), and q = 10 and / = 13 (d). The 
numerically calculated values are given in red circles, the approximated (more accurate) formula is 
represented by blue crosses, and the simple approximated (less accurate) formula is represented by 
a black line. Only odd periods were calculated. Similar results were obtained with up to 12 qubits 
and all possible shifts. 
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Figure 3: (Color online) The average absolute value of the change AG(V') in the entanglement, 
defined in Eq. (|33p . induced by the QFT for a sample of random quantum states (a) and for 
periodic states (b), as a function of the number of qubits, q. The random states are taken from 
a uniform distribution on the 2(7-dimensional complex unit sphere. For the periodic states, the 
average is over all periodic states with a given number of qubits. For the random states the error 
bars represent one standard deviation in each direction, while for the periodic states the distribution 
is extremely narrow, much narrower than the width of the line. 
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